Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome.

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

Biallelic variants in NSUN6 cause an autosomal recessive neurodevelopmental disorder.

PURPOSE: 5-methylcytosine RNA modifications are driven by NSUN methyltransferases. Although variants in NSUN2 and NSUN3 were associated with neurodevelopmental diseases, the physiological role of NSUN6 modifications on transfer RNAs and messenger RNAs remained elusive. METHODS: We combined exome sequencing of consanguineous families with functional characterization to identify a new neurodevelopmental disorder gene. RESULTS: We identified 3 unrelated consanguineous families with deleterious homozygous variants in NSUN6. Two of these variants are predicted to be loss-of-function. One maps to the first exon and is predicted to lead to the absence of NSUN6 via nonsense-mediated decay, whereas we showed that the other maps to the last exon and encodes a protein that does not fold correctly. Likewise, we demonstrated that the missense variant identified in the third family has lost its enzymatic activity and is unable to bind the methyl donor S-adenosyl-L-methionine. The affected individuals present with developmental delay, intellectual disability, motor delay, and behavioral anomalies. Homozygous ablation of the NSUN6 ortholog in Drosophila led to locomotion and learning impairment. CONCLUSION: Our data provide evidence that biallelic pathogenic variants in NSUN6 cause one form of autosomal recessive intellectual disability, establishing another link between RNA modification and cognition.

Bi-allelic TTI1 variants cause an autosomal-recessive neurodevelopmental disorder with microcephaly.

Telomere maintenance 2 (TELO2), Tel2 interacting protein 2 (TTI2), and Tel2 interacting protein 1 (TTI1) are the three components of the conserved Triple T (TTT) complex that modulates activity of phosphatidylinositol 3-kinase-related protein kinases (PIKKs), including mTOR, ATM, and ATR, by regulating the assembly of mTOR complex 1 (mTORC1). The TTT complex is essential for the expression, maturation, and stability of ATM and ATR in response to DNA damage. TELO2- and TTI2-related bi-allelic autosomal-recessive (AR) encephalopathies have been described in individuals with moderate to severe intellectual disability (ID), short stature, postnatal microcephaly, and a movement disorder (in the case of variants within TELO2). We present clinical, genomic, and functional data from 11 individuals in 9 unrelated families with bi-allelic variants in TTI1. All present with ID, and most with microcephaly, short stature, and a movement disorder. Functional studies performed in HEK293T cell lines and fibroblasts and lymphoblastoid cells derived from 4 unrelated individuals showed impairment of the TTT complex and of mTOR pathway activity which is improved by treatment with Rapamycin. Our data delineate a TTI1-related neurodevelopmental disorder and expand the group of disorders related to the TTT complex.

Delineating the Spectrum of Genetic Variants Associated with Bardet-Biedl Syndrome in Consanguineous Pakistani Pedigrees.

This study aimed to find the molecular basis of Bardet-Biedl syndrome (BBS) in Pakistani consanguineous families. A total of 12 affected families were enrolled. Clinical investigations were performed to access the BBS-associated phenotypes. Whole exome sequencing was conducted on one affected individual from each family. The computational functional analysis predicted the variants' pathogenic effects and modeled the mutated proteins. Whole-exome sequencing revealed 9 pathogenic variants in six genes associated with BBS in 12 families. The BBS6/MKS was the most common BBS causative gene identified in five families (5/12, 41.6%), with one novel (c.1226G>A, p.Gly409Glu) and two reported variants. c.774G>A, Thr259LeuTer21 was the most frequent BBS6/MMKS allele in three families 3/5 (60%). Two variants, c.223C>T, p.Arg75Ter and a novel, c. 252delA, p.Lys85STer39 were detected in the BBS9 gene. A novel 8bp deletion c.387_394delAAATAAAA, p. Asn130GlyfsTer3 was found in BBS3 gene. Three known variants were detected in the BBS1, BBS2, and BBS7 genes. Identification of novel likely pathogenic variants in three genes reaffirms the allelic and genetic heterogeneity of BBS in Pakistani patients. The clinical differences among patients carrying the same pathogenic variant may be due to other factors influencing the phenotype, including variants in other modifier genes.

Deleterious, protein-altering variants in the transcriptional coregulator ZMYM3 in 27 individuals with a neurodevelopmental delay phenotype.

Neurodevelopmental disorders (NDDs) result from highly penetrant variation in hundreds of different genes, some of which have not yet been identified. Using the MatchMaker Exchange, we assembled a cohort of 27 individuals with rare, protein-altering variation in the transcriptional coregulator ZMYM3, located on the X chromosome. Most (n = 24) individuals were males, 17 of which have a maternally inherited variant; six individuals (4 male, 2 female) harbor de novo variants. Overlapping features included developmental delay, intellectual disability, behavioral abnormalities, and a specific facial gestalt in a subset of males. Variants in almost all individuals (n = 26) are missense, including six that recurrently affect two residues. Four unrelated probands were identified with inherited variation affecting Arg441, a site at which variation has been previously seen in NDD-affected siblings, and two individuals have de novo variation resulting in p.Arg1294Cys (c.3880C>T). All variants affect evolutionarily conserved sites, and most are predicted to damage protein structure or function. ZMYM3 is relatively intolerant to variation in the general population, is widely expressed across human tissues, and encodes a component of the KDM1A-RCOR1 chromatin-modifying complex. ChIP-seq experiments on one variant, p.Arg1274Trp, indicate dramatically reduced genomic occupancy, supporting a hypomorphic effect. While we are unable to perform statistical evaluations to definitively support a causative role for variation in ZMYM3, the totality of the evidence, including 27 affected individuals, recurrent variation at two codons, overlapping phenotypic features, protein-modeling data, evolutionary constraint, and experimentally confirmed functional effects strongly support ZMYM3 as an NDD-associated gene.

Functional characteristics of a broad spectrum of TBX6 variants in Mayer-Rokitansky-Küster-Hauser syndrome.

PURPOSE: Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) is characterized by congenital absence of the uterus, cervix, and the upper part of the vagina in females. Whole-gene deletion and loss-of-function variants in TBX6 have been identified in association with MRKHS. We aimed to expand the spectrum of TBX6 variants in MRKHS and explore the biological effect of the variant alleles. METHODS: Rare variants in TBX6 were called from a combined multiethnic cohort of 622 probands with MRKHS who underwent exome sequencing or genome sequencing. Multiple in vitro functional experiments were performed, including messenger RNA analysis, western blotting, transcriptional activity assay, and immunofluorescence staining. RESULTS: We identified 16 rare variants in TBX6 from the combined cohort, including 1 protein-truncating variant reported in our previous study and 15 variants with unknown effects. By comparing the prevalence of TBX6 variants in the Chinese MRKHS cohort vs 1038 female controls, we observed a significant mutational burden of TBX6 in affected individuals (P = .0004, odds ratio = 5.25), suggesting a causal role of TBX6 variants in MRKHS. Of the 15 variants with uncertain effects, 7 were shown to induce a loss-of-function effect through various mechanisms. The c.423G>A (p.Leu141=) and c.839+5G>A variants impaired the normal splicing of TBX6 messenger RNA, c.422T>C (p.Leu141Pro) and c.745G>A (p.Val249Met) led to decreased protein expression, c.10C>T (p.Pro4Ser) and c.400G>A (p.Glu134Lys) resulted in perturbed transcriptional activity, and c.356G>A (p.Arg119His) caused protein mislocalization. We observed incomplete penetrance and variable expressivity in families carrying deleterious variants, which indicates a more complex genetic mechanism than classical Mendelian inheritance. CONCLUSION: Our study expands the mutational spectrum of TBX6 in MRKHS and delineates the molecular pathogenesis of TBX6 variants, supporting the association between deleterious variants in TBX6 and MRKHS.

GABBR1 monoallelic de novo variants linked to neurodevelopmental delay and epilepsy.

GABAB receptors are obligatory heterodimers responsible for prolonged neuronal inhibition in the central nervous system. The two receptor subunits are encoded by GABBR1 and GABBR2. Variants in GABBR2 have been associated with a Rett-like phenotype (MIM: 617903), epileptic encephalopathy (MIM: 617904), and milder forms of developmental delay with absence epilepsy. To date, however, no phenotypes associated with pathogenic variants of GABBR1 have been established. Through GeneMatcher, we have ascertained four individuals who each have a monoallelic GABBR1 de novo non-synonymous variant; these individuals exhibit motor and/or language delay, ranging from mild to severe, and in one case, epilepsy. Further phenotypic features include varying degrees of intellectual disability, learning difficulties, autism, ADHD, ODD, sleep disorders, and muscular hypotonia. We functionally characterized the four de novo GABBR1 variants, p.Glu368Asp, p.Ala397Val, p.Ala535Thr, and p.Gly673Asp, in transfected HEK293 cells. GABA fails to efficiently activate the variant receptors, most likely leading to an increase in the excitation/inhibition balance in the central nervous system. Variant p.Gly673Asp in transmembrane domain 3 (TMD3) renders the receptor completely inactive, consistent with failure of the receptor to reach the cell surface. p.Glu368Asp is located near the orthosteric binding site and reduces GABA potency and efficacy at the receptor. GABA exhibits normal potency but decreased efficacy at the p.Ala397Val and p.Ala535Thr variants. Functional characterization of GABBR1-related variants provides a rationale for understanding the severity of disease phenotypes and points to possible therapeutic strategies.

Heterozygous variants in CTR9, which encodes a major component of the PAF1 complex, are associated with a neurodevelopmental disorder.

PURPOSE: CTR9 is a subunit of the PAF1 complex (PAF1C) that plays a crucial role in transcription regulation by binding CTR9 to RNA polymerase II. It is involved in transcription-coupled histone modification through promoting H3K4 and H3K36 methylation. We describe the clinical and molecular studies in 13 probands, harboring likely pathogenic CTR9 missense variants, collected through GeneMatcher. METHODS: Exome sequencing was performed in all individuals. CTR9 variants were assessed through 3-dimensional modeling of the activated human transcription complex Pol II-DSIF-PAF-SPT6 and the PAF1/CTR9 complex. H3K4/H3K36 methylation analysis, mitophagy assessment based on tetramethylrhodamine ethyl ester perchlorate immunofluorescence, and RNA-sequencing in skin fibroblasts from 4 patients was performed. RESULTS: Common clinical findings were variable degrees of intellectual disability, hypotonia, joint hyperlaxity, speech delay, coordination problems, tremor, and autism spectrum disorder. Mild dysmorphism and cardiac anomalies were less frequent. For 11 CTR9 variants, de novo occurrence was shown. Three-dimensional modeling predicted a likely disruptive effect of the variants on local CTR9 structure and protein interaction. Additional studies in fibroblasts did not unveil the downstream functional consequences of the identified variants. CONCLUSION: We describe a neurodevelopmental disorder caused by (mainly) de novo variants in CTR9, likely affecting PAF1C function.

Short arms of human acrocentric chromosomes and the completion of the human genome sequence.

The complete, ungapped sequence of the short arms of human acrocentric chromosomes (SAACs) is still unknown almost 20 years after the near completion of the Human Genome Project. Yet these short arms of Chromosomes 13, 14, 15, 21, and 22 contain the ribosomal DNA (rDNA) genes, which are of paramount importance for human biology. The sequences of SAACs show an extensive variation in the copy number of the various repetitive elements, the full extent of which is currently unknown. In addition, the full spectrum of repeated sequences, their organization, and the low copy number functional elements are also unknown. The Telomere-to-Telomere (T2T) Project using mainly long-read sequence technology has recently completed the assembly of the genome from a hydatidiform mole, CHM13, and has thus established a baseline reference for further studies on the organization, variation, functional annotation, and impact in human disorders of all the previously unknown genomic segments, including the SAACs. The publication of the initial results of the T2T Project will update and improve the reference genome for a better understanding of the evolution and function of the human genome.

The complete sequence of a human genome.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion-base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

CoverageMaster: comprehensive CNV detection and visualization from NGS short reads for genetic medicine applications.

CoverageMaster (CoM) is a copy number variation (CNV) calling algorithm based on depth-of-coverage maps designed to detect CNVs of any size in exome [whole exome sequencing (WES)] and genome [whole genome sequencing (WGS)] data. The core of the algorithm is the compression of sequencing coverage data in a multiscale Wavelet space and the analysis through an iterative Hidden Markov Model. CoM processes WES and WGS data at nucleotide scale resolution and accurately detects and visualizes full size range CNVs, including single or partial exon deletions and duplications. The results obtained with this approach support the possibility for coverage-based CNV callers to replace probe-based methods such as array comparative genomic hybridization and multiplex ligation-dependent probe amplification in the near future.

Biallelic truncation variants in ATP9A are associated with a novel autosomal recessive neurodevelopmental disorder.

Intellectual disability (ID) is a highly heterogeneous disorder with hundreds of associated genes. Despite progress in the identification of the genetic causes of ID following the introduction of high-throughput sequencing, about half of affected individuals still remain without a molecular diagnosis. Consanguineous families with affected individuals provide a unique opportunity to identify novel recessive causative genes. In this report, we describe a novel autosomal recessive neurodevelopmental disorder. We identified two consanguineous families with homozygous variants predicted to alter the splicing of ATP9A which encodes a transmembrane lipid flippase of the class II P4-ATPases. The three individuals homozygous for these putatively truncating variants presented with severe ID, motor and speech impairment, and behavioral anomalies. Consistent with a causative role of ATP9A in these patients, a previously described Atp9a-/- mouse model showed behavioral changes.

Three decades of the Human Genome Organization.

The Human Genome Organization (HUGO) was initially established in 1988 to help integrate international scientific genomic activity and to accelerate the diffusion of knowledge from the efforts of the human genome project. Its founding President was Victor McKusick. During the late 1980s and 1990s, HUGO organized lively gene mapping meetings to accurately place genes on the genome as chromosomes were being sequenced. With the completion of the Human Genome Project, HUGO went through some transitions and self-reflection. In 2020, HUGO (which hosts a large annual scientific meeting and comprises the renowned HUGO Gene Nomenclature Committee [HGNC], responsible for naming genes, and an outstanding Ethics Committee) was merged with the Human Genome Variation Society (HGVS; which defines the correct nomenclature for variation description) and the Human Variome Project (HVP; championed by the late Richard Cotton) into a single organization that is committed to assembling human genomic variation from all over the world. This consolidated effort, under a new Executive Board and seven focused committees, will facilitate efficient and effective communication and action to bring the benefits of increasing knowledge of genome diversity and biology to people all over the world.

Immune Dysregulation and the Increased Risk of Complications and Mortality Following Respiratory Tract Infections in Adults With Down Syndrome.

The risk of severe outcomes following respiratory tract infections is significantly increased in individuals over 60 years, especially in those with chronic medical conditions, i.e., hypertension, diabetes, cardiovascular disease, dementia, chronic respiratory disease, and cancer. Down Syndrome (DS), the most prevalent intellectual disability, is caused by trisomy-21 in ~1:750 live births worldwide. Over the past few decades, a substantial body of evidence has accumulated, pointing at the occurrence of alterations, impairments, and subsequently dysfunction of the various components of the immune system in individuals with DS. This associates with increased vulnerability to respiratory tract infections in this population, such as the influenza virus, respiratory syncytial virus, SARS-CoV-2 (COVID-19), and bacterial pneumonias. To emphasize this link, here we comprehensively review the immunobiology of DS and its contribution to higher susceptibility to severe illness and mortality from respiratory tract infections.

History of the methodology of disease gene identification.

The past 45 years have witnessed a triumph in the discovery of genes and genetic variation that cause Mendelian disorders due to high impact variants. Important discoveries and organized projects have provided the necessary tools and infrastructure for the identification of gene defects leading to thousands of monogenic phenotypes. This endeavor can be divided in three phases in which different laboratory strategies were employed for the discovery of disease-related genes: (i) the biochemical phase, (ii) the genetic linkage followed by positional cloning phase, and (iii) the sequence identification phase. However, much more work is needed to identify all the high impact genomic variation that substantially contributes to the phenotypic variation.

Phenotypic expansion of CACNA1C-associated disorders to include isolated neurological manifestations.

PURPOSE: CACNA1C encodes the alpha-1-subunit of a voltage-dependent L-type calcium channel expressed in human heart and brain. Heterozygous variants in CACNA1C have previously been reported in association with Timothy syndrome and long QT syndrome. Several case reports have suggested that CACNA1C variation may also be associated with a primarily neurological phenotype. METHODS: We describe 25 individuals from 22 families with heterozygous variants in CACNA1C, who present with predominantly neurological manifestations. RESULTS: Fourteen individuals have de novo, nontruncating variants and present variably with developmental delays, intellectual disability, autism, hypotonia, ataxia, and epilepsy. Functional studies of a subgroup of missense variants via patch clamp experiments demonstrated differential effects on channel function in vitro, including loss of function (p.Leu1408Val), neutral effect (p.Leu614Arg), and gain of function (p.Leu657Phe, p.Leu614Pro). The remaining 11 individuals from eight families have truncating variants in CACNA1C. The majority of these individuals have expressive language deficits, and half have autism. CONCLUSION: We expand the phenotype associated with CACNA1C variants to include neurodevelopmental abnormalities and epilepsy, in the absence of classic features of Timothy syndrome or long QT syndrome.

PIGG variant pathogenicity assessment reveals characteristic features within 19 families.

PURPOSE: Phosphatidylinositol Glycan Anchor Biosynthesis, class G (PIGG) is an ethanolamine phosphate transferase catalyzing the modification of glycosylphosphatidylinositol (GPI). GPI serves as an anchor on the cell membrane for surface proteins called GPI-anchored proteins (GPI-APs). Pathogenic variants in genes involved in the biosynthesis of GPI cause inherited GPI deficiency (IGD), which still needs to be further characterized. METHODS: We describe 22 individuals from 19 unrelated families with biallelic variants in PIGG. We analyzed GPI-AP surface levels on granulocytes and fibroblasts for three and two individuals, respectively. We demonstrated enzymatic activity defects for PIGG variants in vitro in a PIGG/PIGO double knockout system. RESULTS: Phenotypic analysis of reported individuals reveals shared PIGG deficiency-associated features. All tested GPI-APs were unchanged on granulocytes whereas CD73 level in fibroblasts was decreased. In addition to classic IGD symptoms such as hypotonia, intellectual disability/developmental delay (ID/DD), and seizures, individuals with PIGG variants of null or severely decreased activity showed cerebellar atrophy, various neurological manifestations, and mitochondrial dysfunction, a feature increasingly recognized in IGDs. Individuals with mildly decreased activity showed autism spectrum disorder. CONCLUSION: This in vitro system is a useful method to validate the pathogenicity of variants in PIGG and to study PIGG physiological functions.